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摘 要 : 针对 传统 哈里 斯 鹰 优化 算法 (Harris Hawk Optimization Algorithm，HHO) 求 解 精 度 低 、 迭 代 后 期 收敛 慢 以 及 
易 早 熟 收 化 的 不 足 ， 提 出 一 种 融合 随机 无 迹 0 变异 的 改进 哈里 斯 座 优 化 算法 (Opposite-Support Harris Hawk 
Optimization Algorithm，OSHHO)。 利 用 伪 对 立 和 伪 反 射 学 习 机 制 加 快 HHO 算法 的 寻 优 效率 ， 提 升 全 局 收敛 能 
提出 能 量 因子 非 线性 周期 性 调整 机 制 均 衡 算 法 的 全 局 搜索 与 局 部 开发 能 力 ; 设计 随机 无 迹 变异 机 制 使 算法 能 够 跳 离 
局 部 极 值 点 ， 进 一 步 提 高 算法 的 寻 优 精度 和 收敛 速度 。 在 八 个 基准 函数 上 测试 了 OSHHO 工法 的 寻 优 性 能 。 为 了 验 
证 算法 的 实用 性 ， 利 用 改进 哈里 斯 应 优化 算法 OSHHO 优化 支持 向 量 机 模型 (Support Vector Machine，SVM) 的 关键 
参数 ， 并 将 模型 应 用 于 新 冠 肺炎 疫情 传播 预测 问题 。 随 机 选取 若干 天 新 冠 肺炎 疫情 日 增 确认 人 数 作为 样本 数据 ， 利 
用 训练 好 的 SVM 模型 预测 新 冠 肺 类 疫情 变化 趋势 。 结 果 证 实 算 法 具有 更 高 的 预测 精度 和 更 快 的 收敛 速度 。 
关键 词 : 哈里 斯 优化 算法 ; 伪 对 立 学 习 ; 伪 反射 学 习 ; 无 迹 变异 ; 传播 预测 ; 支持 向 量 机 ; 新 冠 肺炎 
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Improved Harris hawk optimization algorithm based on random traceless o mutation and ITS applications 
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Mathematics & Informatics/College of Software, South China Agricultural University, Guangzhou 510642, China) 


, Abstract: Aiming at the shortcomings of the traditional harris hawk optimization algorithm, such as low precision, slow 
LO convergence speed in the later iteration and fall into easily a premature convergence. An improved harris hawk optimization 
algorithm OSHHO integrating with random traceless mutation is proposed. The pseudo back and pseudo reflection learning 
mechanism is used to speed up the optimization efficiency of HHO algorithm and improve the global convergence ability. The 
nonlinear periodic adjustment mechanism of energy factor ls proposed to improve the global search and local development 
ability of the equalization algorithm. The random traceless point mutation mechanism is designed to make the algorithm jump 


和 away from the local extreme points, and further improve the optimization accuracy and convergence Speed of the algorithm . 
Eight benchmark functions are used to test the optimization performance of OSHHO algorithm. And two typical engineering 
design problems are solved by OSHHO to further versify the feasibility and the efficiency of the proposed algorithm. The 


results show that the new algorithm has higher prediction accuracy and faster convergence Speed. 
Key words: Harris hawk algorithm; pseduo opposition-learning; pseduo reflection-learning; traceless mutation; diffuse 


prediction; support vector machine; covid-19 

0 引言 与 局 部 开发 过 程 的 不 足 。 为 此 ， 出 现 了 许多 针对 标准 哈里 斯 
攻 鹰 优化 算法 的 改进 算法 ， 这 类 算法 通常 针对 某 一 方面 来 增强 
单 体 智 能 优化 算法 是 通过 模拟 自然 界 物 种 的 社会 行为 机 其 局 部 开发 能 力 或 全 局 搜索 能 力 。 一 种 思路 是 融合 其 他 智能 
制 衍生 出 的 一 种 优化 方法 ， 这 类 算法 具有 适应 性 、 自 学 习 以 ” 算法 对 其 进行 改进 ， 如 文献 [10] 将 差分 进化 算法 DE 融入 
及 群体 搜索 性 特征 ， 在 解决 高 维 复杂 和 多 极 值 优化 问题 时 ， HHO 算法 ， 将 种 群 分 割 为 两 个 部 分 分 别 以 HHO 和 DE 算法 
拥有 传统 方法 不 可 比拟 的 优势 。 近 年 来 ， 已 经 涌现 出 鲸鱼 优 ”进行 个 体 进化 ， 能 够 提高 寻 优 精度 。 文 献 [1H] 将 正 余 弦 算 法 
化 算法 (Whale Optimization Algorithm,WOA)JI、 灰 狼 优 化 算 ”SCA 融入 HHO 算法 ， 设 计 了 混合 变异 算法 ， 一 定 程 度 避 免 
法 (Grey Wolf Optimization,GWO) 中 、 蝴 蝶 优 化 算法 (Butterfly ”了 HHO 算法 陷入 局 部 最 优 ， 算 法 搜索 能 力 得 到 了 改进 。 文 
Optimization Algorithm,BOA)B1 和 遵 海 鞘 算 法 (Salp Swarm [12] 则 将 模拟 退火 算法 SSA 和 HHO 融合 ， 以 此 提高 HHO 
Algorithm,SSA) 负 等 一 批 寻 优 能 力 较 强 的 群体 智能 算法 。 而 哈 。 算法 的 全 局 收敛 能 力 ， 并 以 参数 寻 优 问题 进行 了 验证 。 男 一 种 
斯 雇 优 化 算法 (Harris Hawks Optimization,HHOJDG] 是 在 2019 ”思路 是 利用 相关 方法 对 HHO 的 结构 进行 改进 ， 如 文献 [13] 在 
学 者 Heidari 等 提出 的 一 种 群体 智能 优化 算法 ， 模 拟 了 ”HHO 算法 的 搜索 阶段 利用 柯 西 变异 机 制 提高 种 群 多 样 性 ， 并 利 

自然 界 中 哈里 斯 鹰 的 群体 捕食 行为 ， 以 种 群 对 猪 物 的 包围 、 随机 收缩 改进 能 量 方程 ， 均 衡 全 局 搜索 和 局 部 开发 过 程 。 文 
攻击 过 程 建 立 算法 的 数学 模型 。HHO 算法 综合 性 能 较 优 、 结 献 [14] 引 入 方形 邻 域 拓扑 结构 改进 种 群 个 体 组 成 ， 通 过 纵横 双 
构 简 单 ， 同 时 依赖 参数 少 ， 一 经 提出 ， 已 经 广泛 应 用 于 诸多 ] 随 机 况 食 提高 HHO 算法 的 开发 能 力 。 文 献 [15] 引 入 伪 对 立 学 
领域 , 如 : 图 像 处 理 器 、 神 经 网 络 训练 四、 电力 控制 名、TDOA ” 习 提 高 HHO 算法 初始 种 群 的 多 样 性 ， 进 一 步 提升 了 算法 的 收 
定位 外 等 ,然而 ,标准 HHO 算法 本 身 还 是 存在 着 易 早 熟 收敛 、 ”人 钱 速 度 。 文献 [16] 设 计 基于 对 立 学 习 和 正 交 学 习 的 改进 HHO 
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目前 改进 HHO 算法 主要 针 
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rs 为 (0,1) 间 的 随机 量 。 
b) 硬 包围 。 若 4>0.5 且 |Bl<0.5 时 ， 表 明 猫 物 没有 足够 能 

庶 会 迅速 捕食 猎物 。 其 数学 模型 为 

X(t+1)= Xt) —E|AX() (5) 
c) 渐 近 快 速 俯冲 式 软 包 围 。 若 4<0.5 且 |El20.5 时 ， 表 明 

猎物 仍 有 能 量 逃 脱 ， 但 哈里 斯 雇 种 群 会 环绕 、 包 围 并 追捕 猪 

物 ， 该 数学 模型 为 

Y: X00) —E|IX,00(t)) — XD FY) <FX,O)) (0) 

Z:Y+SxLF(D), F(Z)<F(X,()) 


wll 
区 
上 
|T 阅 


各 


xorD-| 


的 关 


键 参数 进行 寻 优 ， 建 立 一 种 面向 新 冠 


哈里 斯 鹰 优 化 算法 HHO 


斯 鹰 优化 算法 HHO 是 一 种 通 


协作 捕食 行为 而 建立 的 随机 式 群 


宗 、 围 捕 、 攻 击 对 猎物 进行 


索 至 了 


鹰 个 体 i 的 原 位 
迭代 + 时 种 群 中 的 全 局 
随机 量 ，[15,ub] 为 个 
代 t 时 种 群 的 平均 位 


其 中 : NN 为 哈日 
1.2 搜索 至 开发 的 过 渡 切 换 

HHO 算法 实现 从 全 
逃逸 能 量 因子 控制 


这 一 阶段 中 ， 


= 折 谁 个 体 (位置 ) 可 
Pp 逐步 靠近 猎物 ， 即 最 优 解 


视 为 问 
。 算 法 可 以 分 为 


ff 发 间 的 切换 以 及 所 
ee 


EE ，Xi(t+1) 为 个 


采 三 个 阶段 。 


本 i 的 更 新 位 


过 观察 哈 
利索 算法 。 该 种 群 可 以 通过 
高 效 协 作 捕食 。HHO 算法 中 ， 
题 的 一 个 候选 解 ， 捕 食 过 
全 局 搜索 、 搜 


， 克 (办 为 哈里 斯 
二 和 oo0t( 人 为 


最 优 个 体 ，g、ri(i=1,2,3,4) 为 (0,1) 间 的 


里 斯 座 群 的 


斯 座 会 以 两 种 等 概率 的 随机 行为 方式 
有 索 猫 物 ， 休学 大 


trD=txwO-X OO -abrnb -bg<0s (1) 


其 中 : 闷 oa(0) 为 和 欠 代 上 时 随机 选择 个 体 的 位 


其 中 : 为 问题 维度 ,5 为 大 小 为 Dxl 的 (0,1) 间 随机 行 向 量 ， 
LF() 为 服从 Levy 飞行 分 布 的 函数 ， 计 算 方法 为 
HAxI j T(U+B)xsin(zB/2) 
LF(x)=0.01 ne (Tarp DBI) (7) 
其 中 : 4、Y 为 (0,1) 间 随机 量 ，p=1.5。 

d) 渐 近 快 速 俯冲 式 硬 包围 。 若 4<0.5 且 |B<0.5 时 ， 表 明 
猪 物 没有 充足 能 量 逃 逸 ， 哈 里 斯 雇 种 群 会 环绕 建立 一 个 硬 包 
捕食 猎物 ， 该 数学 模型 为 
Y :Xam (D)—E|TXn) —X, (DFOY) < F(X) 
Z:Y+SxLF(D), F(Z)<F(X.,()) (8) 


xcrD-| 


2 ”改进 哈里 斯 雇 优 化 算法 OSHHO 


2.1 伪 对 立 和 伪 反 射 学 习 机 制 

对 立 学 习 (Opposition-based Learning,OBL) 可 以 通过 同步 
考虑 候选 个 体 及 其 对 立 解 的 方式 ， 通 过 择优 保留 至 下 一 代 种 
群 的 方式 提升 种 群 多 样 性 ， 从 而 提升 智能 优化 算法 的 求解 精 
度 和 收敛 速度 。 相 关 模 型 为 : 令 d 维 空间 内 个 体 i 的 位 置 
了 一 (XiIyXi2.…Xiq)，XijS[1bsubj]，[1bjsubj] 为 个 体 在 j 维 空 间 取 
值 范 围 ， 庆 1,2,…,d。 个 体 i 的 对 立 点 为 XY 守 (xil ,Xi2.… id )， 
个 体 i 的 伪 对 立 点 为 天”=(xil1 ”xi2 ”xia”), 个 体 i 的 伪 反 射 
点 为 名 ”=CYil ”Xi2”.….Xid”)， 定 义 分 别 为 


本 搜索 空间 下 界 与 


其 中 : Eo 为 猎物 初始 的 逃逸 
为 当前 迁 代 次 数 ,Ts 为 最 大 迭代 次 数 。 切换 原理 
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段 中 , HHO 算法 可 以 按 四 


搜索 ， 包 括 软 包围 、 


. 围 、 渐 近 快 速 俯 六 


快速 俯冲 式 硬 包 围 。 
为 猎物 成 功 逃 脱 ; 若 4>0.5， 则 猎物 逃脱 失败 。 和 逃逸 能 


五 表示 哈里 
若 |E|<0.5， 则 执行 硬 包围 
a) 软 包围 。 若 120.5 目 


令 1 为 猎物 的 逃 


斯 鹰 的 软 硬 攻 击 策略 , 若 


局 部 开采 。 


种 不 同 策略 对 猎物 进行 
PF 式 软 包围 和 渐 近 
脱 概率 ， 若 4<0.5， 则 认 


如 >0.5， 则 执行 软 
此 共有 以 下 四 种 捕食 策 


为 : 阁 |E1， 


大 二 


>0.5 时 ， 表 明 猫 物 拥有 足够 能 

量 逃 脱 ， 哈 里 斯 座 会 环绕 包围 猎物 ， 逐 渐 消 耗 猎物 能 量 ， 在 
物 失 去 能 量 时 完成 猎物 捕食 。 其 数学 模型 为 

1 —E|JX,m 0) — XD| @) 


AX (1) = Xappi 


其 中 : 人 入 XD) 为 猎物 与 哈 


一 Xi(D) 


里 斯 鹰 个 体位 置 的 间距 ，, 户 2(1-75)， 


Xj =Ib; +ubj; — x (9) 
Xj "=rand[(lb; +ub))/2,x, ] (10) 
Xj" =rand[x;;, (lb; +ub,)/2] (1 1) 


在 二 维 空间 中 ， 个 体 下 (x1,x2) 的 伪 对 立 点 和 伪 反 射 点 的 
取 值 范围 如 图 1 所 示 。 伪 反 射 点 X=(x1”,x2”) 的 取 值 范围 为 
区 域 4, 伪 对 立 点 子 ”=(x1”,x2”) 的 取 值 范围 为 区 域 B。 根据 图 
示 结 果 ， 伪 反射 点 邓 ” 比 伪 对 立 点 卫 ”" 更 接近 于 候选 解 ， 可 以 
在 候选 解 的 邻近 区 域 进行 充分 的 局 部 开发 。 伪 对 立 点 离 候选 
解 位 置 较 远 ， 则 可 以 实现 更 广泛 的 全 局 搜索 ， 开 辟 候 选 解 未 
探索 的 空间 。 
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图 1 伪 对 立 点 和 伪 反 射 点 在 搜索 空间 内 的 分 布 
Fig.1 Distribution ofpseudo opposition points and pseudo reflection 
points in Search Space 

对 于 HHO 算法 而 言 ， 当 | 如 >0.5 且 4>0.5 时 , 猎物 拥有 足 
够 能 量 以 跳跃 式 逃 脱 包围 ， 这 使 得 哈里 斯 鹰 个 体 需要 扩大 搜 
索 范围 来 捕获 猎物 。 为 此 ， 在 软 包 围 时 期 的 位 置 更 新 中 引入 
伪 对 立 学 习 机 制 。 根 据 图 1 所 示 结 果 ， 伪 对 立 点 离 候 选 解 的 
位 置 较 远 ,在 远离 候选 解 的 区 域内 生成 当前 解 的 伪 对 立 解 后 ， 
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个 体能 够 实现 更 广泛 的 全 局 搜索 ， 拓 展 候选 解 未 搜索 过 的 区 
域 ， 以 此 增强 该 过 程 中 的 种 群 多 样 性 ， 并 最 后 通过 贪 禁 选择 
策略 保留 原 个 体 和 伪 对 立 解 个 体 进入 下 一 代 种 群 。 令 
天 old(tt1) 为 HHO 算法 软 包 围 时 的 位 置 更 新 ， 则 引入 伪 对 立 


学 习 后 的 位 置 更 新 方式 为 
X44D)= 人 G+D, Xialt +D))< /0 "(t+D)) (12) 
Xi CC+D, f (Xioalt +D))> F(X "t+1)) 
其 中 : Kiva lt +1)=AX,(D)—E|JX,d) — XD| (13) 
Xi "(t+h)= om (+h),..., Kia.oa "(t+1)) (14) 
Kjoa t+) = rand[(lb; +ub))/2,lb; +ub) —xi jaa (t+)] (15) 


当 ||<0.5 且 4>0.5 时 ， 猪 物 没有 足够 能 量 逃 脱 ， 这 使 得 
哈里 斯 鹰 个 体会 在 小 范围 内 迅速 捕食 猎物 。 为 此 ， 在 硬 包围 
时 期 的 位 置 更 新 中 引入 伪 反 射 学 习 机 制 。 根 据 图 1 所 示 结 果 ， 
反射 点 更 加 接近 候选 解 ， 在 邻近 候选 解 的 区 域内 生成 当前 
擅 反 射 解 后 ， 个 体能 够 实现 更 强 的 局 部 开发 能 力 ， 在 
空间 内 实现 精细 搜索 ， 并 通过 贪 焚 选 择 策 略 保留 原 个 体 
和 伪 反 射 解 个 体 进入 下 一 代 种 群 。 引 入 伪 反 射 学 习 后 的 位 


| < f Xi "HD) 6 
XHD, fF Xi (t+) > FX "(+1)) 
其 中 : Kia (t+D) = X00) -EIAX,D| (17) 
Xi "HD = G0 (+ Ds a "(t+D)) (18) 
Kjaa (+1) = rand[x jaa (t+1),(1b; +ub) /2 (19) 


通过 在 局 部 开发 过 程 中 引入 伪 对 立 学 习 和 伪 反 射 学 习 机 
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A 1 
= Coss " | (24) 


max 


中 ,i=1,2,.…,4，SDx 为 协 方差 矩阵 ， re 为 控制 因子 ， 用 于 控 


制 在 最 优 解 邻近 区 域内 的 开采 过 程 。OSHHO 算法 实施 个 体 
变异 扰动 后 ， 将 从 生成 的 变异 个 体 中 选择 最 优 个 体 作为 新 的 
种 群 最 优 解 ， 可 表示 为 

Xm (t+D) =argmin{f(X;)} (25) 

_A 
机 rw 
6 
二 


图 2 d=3 时 6 点 分 布 
Fig.2 Distribution of o points when d=3 
2.3 能 量 因子 非 线性 调整 机 制 
对 于 HHO 这 类 群体 智能 优化 算法 而 言 ， 全 局 搜索 与 局 
部 开发 过 程 的 均衡 切换 对 算法 寻 优 极为 关键 。 算 法 迭代 早期 ， 
需要 更 强大 的 全 局 搜索 能 力 ， 以 便 维 持 种 群 分 布 的 多 样 性 ; 
而 在 迭代 后 期 ， 则 需要 更 好 的 局 部 开发 能 力 ， 以 保证 算法 在 


制 ， 根 据 个 体 不 同 的 搜索 策略 ， 能 够 加 快 HHO 算法 的 寻 优 
效率 ， 提 升 全 局 收敛 能 力 。 
2.2 ”随机 无 迹 点 变异 机 制 

在 HHO 算法 欠 代 寻 优 过 程 中 ， 和 迭代 后 期 种 群 个 体 均 会 
向 着 全 局 最 优 解 各 wwwx 所 在 区 域 逼 近 ， 这 会 导致 种 群 多 样 性 
缺失 ， 和 迭代 停滞 进而 得 到 局 部 最 优 解 ， 这 也 是 智能 优化 算法 
的 固有 不 足 。 为 此 ，OSHHO 算法 引入 一 种 针对 当前 最 优 解 
总 ooiz 的 随机 无 迹 o 点 变异 机 制 , 使 算法 具备 跳 离 局 部 最 优 的 


局部 范围 内 的 精细 开采 ， 加 快 算法 收 化 。 由 HHO 算法 执行 
过 程 可 知 ， 逃 逸 能 量 因子 一 决定 了 HHO 算法 的 全 局 搜索 到 
局 部 开采 之 间 的 切换 ， 是 衡量 HHO 算法 寻 优 能 力 的 重要 指 
标 。 当 | 如 >0.5 时 , 哈里 斯 硬 个 体 进行 软 包 围 ( 渐 近 快速 俯冲 式 
软 包 围 ), 种 群 将 扩大 搜索 范围 , 对 应 于 算法 的 全 局 搜索 阶段 ; 
当 | 如 <0.5 时 , 哈里 斯 硬 个 体 进行 硬 包 围 ( 渐 近 快速 俯冲 式 硬 包 
姑 )， 种 群 将 缩小 搜索 范围 ， 对 应 于 算法 的 局 部 开发 阶段 。 然 
， 式 (3) 表 明 , EE 值 在 迭代 过 程 中 是 线性 单 周 期 递减 的 ， 这 


mn 
att 


能 力 。 无 迹 o 点 变换 主要 通过 随机 变量 均值 xave 及 其 方差 值 
SD 产生 24+1 个 随机 点 ， 用 于 评估 随机 变量 函数 y=g(x) 的 均 


值 和 方差 。 对 称 无 迹 c 点 的 生成 方式 为 
G0 "= Xa 
Gi'=Xwe + (Vd +K)SD.,),i=1,2,...,d (20) 


Gia = Xow ~ (Nd +K)SD.)i,i=1,2,...,d 


其 中 ，x=a?(4+4)-4，4=3-d，a 表示 极 小 正 数 ，4 为 问题 空间 
E 度 ，xave 为 变量 x 的 均值 ，SD: 为 变量 x 的 协 方差 矩阵 ，( 
表示 (d+x)SD: 的 平方 根 矩 阵 的 第 i 列 。 无 迹 o。 点 通过 函数 
y 产 g(&7) 的 传播 ， 可 以 生成 24+1 个 y 值 。 则 随机 变量 y 的 均 
值 yave 和 方差 SD， ee 


2d+l 


Yo 4H (21) 


A 


1 她 
5D, = 2 (1 — Yang NF — ag (22) 


如 图 2 为 三 维 空间 中 , 即 d=3 时 无 迹 o 点 的 分 布 示 意图 。 
利用 随机 无 迹 o 点 对 HHO 算法 中 的 全 局 最 优 解 apvi(?) 


对 其 进行 随机 无 迹 变 换 , 在 最 优 解 和 Zavvi( 四 的 邻近 区 域内 作 进 
一 步 高 质量 局 部 开发 ， 提 升 HHO 算法 的 寻 优 精度 。 对 最 优 


进行 变异 的 主要 思路 是 :将 最 优 解 Xrapbi( 四 所 处 位 置 作为 均值 ， 


并 不 符合 哈里 斯 鹰 种 群 需要 多 轮 次 协同 围捕 猎物 的 自然 规律 ， 

线性 递减 的 5 值 会 使 得 每 次 迭代 中 的 能 量变 化 和信 X(D) 是 
值 。 为 此 ，OSHHO 算法 引入 一 种 针对 能 量 因子 的 非 线 性 
期 性 调整 策略 ， 描 述 哈 里 斯 诬 围 捕猎 物 的 过 程 。 具 体 地 ， 
因子 的 更 新 公式 中 以 对 数 函 数 掺 述 非 线性 周期 性 ， 具 体 
定义 为 


全 珂 斑 杰 


E=2Ell-In(d+ ze —1)3)] (26) 


中 ,1 为 当前 迭代 次 数 ，Tmax 为 最 大 迭 代 次 数 ，e 为 自然 常 
Eo 为 猎物 的 初始 能 量 。 根据 式 (26), 迭代 早期 , E 的 衰减 
速度 较 慢 , 迭代 后 期 , 三 的 衰减 速度 变 快 。 如 此 , 迭代 早期 ， 
OSHHO 算法 的 种 群 能 够 进行 更 加 充分 的 全 局 搜索 ， 最 大 限 
度 地 提升 种 群 多 样 性 ;迭代 后 期 ，OSHHO 算法 能 够 更 快 的 
收敛 ， 从 而 更 均衡 的 实现 欠 代 早期 全 局 搜索 与 迭代 后 期 局 部 
开发 间 的 稳定 切换 , 进一步 提高 算法 的 寻 优 精度 和 收敛 速度 
3 为 OSHHO 算法 的 执行 流程 。 


3 ”基于 OSHHO 算法 优化 支持 向 量 机 的 新 冠 肺炎 
疫情 传播 预测 模型 


3.1 支持 向 量 机 SVM 预测 新 冠 肺炎 疫情 传播 
2019 年 年 末 爆 发 的 新 冠 肺炎 让 中 国 居民 始 料 未 及 , 该 传 


地 装 


解 已 copz( 人 0 进行 c 点 变异 ， 生 成 24+1 个 随机 变异 个 体 ， 变 异 
个 体 如 下 : 

Xo'= Xp (i) 

X= Xt) + (Vd + OSDy), (23) 
Xi "= Xi (D1(Vd + x)SDx ), 


染病 具有 传播 感染 快 、 死 亡 率 高 、 持 续 时 间 长 和 防疫 难度 大 
的 特点 ， 对 国内 经 济 及 医疗 卫生 体系 造成 了 巨大 冲击 。 在 抗 
击 新 冠 肺炎 过 程 中 ， 除 生物 医疗 工作 者 在 加 紧 研 制 疫苗 ， 相 
关 研 究 人 员 也 在 通过 利用 数学 模型 对 疫情 传播 预测 进行 研究 。 
如 : 文献 [18] 针 对 疫情 首 个 爆发 城市 武汉 封 城 前 后 的 疫情 传 
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播 模型 进行 分 析 ， 通 过 差分 递 推 和 SIR 模型 建立 疫情 传播 的 步骤 4: 初始 化 哈里 斯 座 种 群 {Xi 二 1,2,…,N}， 将 个 体位 
预测 模型 。 文 献 [19] 利 用 次 指数 增长 拟 合 疫 情 爆 发 早期 状态 ， 置 编码 为 (c, 0)。 
并 对 传播 模式 进行 建 模 。 文 献 [20] 将 长 短期 记忆 网 络 LSTM 聚 5: 计算 个 体 适应 度 ， 确 定 当前 最 优 解 。 
和 SEIR 模型 融合 ， 建 立 了 针对 新 冠 肺 炎 疫 情 传 播 的 预测 模 步骤 6: 对 当前 最 优 解 进行 无 迹 o 变异 , 并 更 新 最 优 解 。 
型 LS-Net, 其 预测 精度 优 于 传统 SEIR 模型 。 文献 [21] 建 立 了 步骤 7: 根据 式 (26) 更 新 猎物 的 逃逸 能 量 因 子 E。 
基于 极端 梯度 提升 树 XGBoost 的 新 冠 肺炎 疫情 传播 预测 模 步骤 8: 若 |E21 且 g>0.5， 利 用 式 (1) 第 1 个 公式 更 新 个 体 
型 ， 模 型 具有 良好 的 预测 能 力 。 受 此 启发 ， 本 文 的 目标 是 结 位 置 ; 若 |EP1 且 gq<0.5， 利 用 式 (1) 第 2 个 公式 更 新 个 体位 置 。 
合群 体 智能 优化 算法 和 支持 向 量 机 SVM 模型 ， 建 立新 冠 肺 步骤 9: 若 |&|20.5 且 4>0.5， 先 利用 式 (4) 更 新 个 体位 置 ， 
炎 疫 情 传播 的 预测 模型 。 并 利用 式 (10) 计 算 其 伪 对 立 点 位 置 , 再 根据 式 (12) 确 定 个 体位 置 
CO 步骤 10: 若 |lEl20.5 且 4<0.5， 利 用 式 (6) 更 新 个 体位 置 。 
i 步骤 11: 若 |2I<0.5 且 人 2>0.5， 先 利用 式 (5) 更 新 个 体位 置 ， 
| 并 利用 式 (11) 计 算 其 伪 反 射 点 位 置 , 再 根据 式 (16) 确 定 个 体位 置 。 
计算 种 群 个 体 适应 度 ， 确 定 最 优 解 步骤 js 若 |E <0.5 且 1<0.5， 利 ) j 式 (8) 更 新 个 体位 置 
ET Ks 13: 更 新 全 局 最 优 解 及 其 适应 度 值 ; 判断 最 优 个 体 
位 置 是 否 超过 (c， 中 的 边界 上 限 ， 若 出 现 越界 ， 则 以 相应 位 
sg 置 上 下 限 修正 个 体位 
， < 步 对 14， 判断 算 法 终止 条 件 ， 若 满足 ， 返 回 至 步 曼 7 执 
| 后 5 到 行 ; 否则 ， 停 止 算 法 闪 代 ， 输 出 最 优 解 Y=(c*,o”)。 
一 < i 步骤 15: 以 了 Y=(c*,o”) 对 SVM 模型 初始 化 ， 并 利用 训练 
人 证 | | | ] | 样本 和 测试 样本 检测 预测 精度 ， 得 到 新 冠 肺炎 确诊 预测 值 。 
| 图 4 是 预测 模型 的 流程 图 。 
| I I I 个 体 编码 ” 怀 一 一 和 
一 下 终归 种 群 初始 化 人 es 和 可 本 委 和 
| 及 参 
人 计算 适应 度 一 一 一 > 得 到 参数 最 优 解 
下 确定 最 优 解 计算 误差 
图 3 OSHHO 算法 的 执行 流程 根 掺 GOSHHO 算 沃 的 子 | | 
Fig. 3 Flow of OSHHO 代 机 制 更 产 个 体位 轩 To 
根据 SVM 模型 可 知 ， 惩 罚 参数 c 和 宽度 参数 o 是 影响 预 CE 
测 模 型 最 关键 的 两 个 参数 。 参 数 c 表示 预测 模型 复杂 度 与 误差 是 
逼近 的 折 中 因子 ， 参 数 o 用 于 控制 核 函 数 的 径 向 范围 。 若 c 值 输出 最 优 解 输出 新 冠 肺炎 确诊 预测 结果 
较 大 ， 则 便于 实现 更 高 的 数据 拟 合 度 ， 但 相应 会 降低 模型 的 泛 图 4 OSHHO 算法 优化 SVM 的 新 冠 肺炎 疫情 传播 预测 模型 
化 能 力 ， 反 之 亦 然 ， 而 若 o 值 较 大 ， 则 会 降低 模型 的 复杂 度 和 Fig. 4 Covid-19 SVM prediction model for pneumonia epidemic 
数据 的 线性 可 分 割 度 ， 较 小 的 o 值 则 可 能 导致 数据 的 严重 过 拟 situation based on OSHHO algorithm 
合 ， 预 测 模型 复杂 度 也 将 趋 近 无 穷 。 由 于 两 个 参数 对 模型 的 泛 4 ”实验 分 析 
化 能 力 和 预测 精度 都 有 着 直接 影响 ， 新 冠 疫情 预测 的 SVM 模 Be 
型 无 论 是 复杂 度 ， 还 是 其 泛 化 能 力 都 取决 于 能 和 否 求解 最 优 的 (c， 4.1 实验 1: 针对 OSHHO 算法 的 数值 仿真 
四 对 。 标准 SVM 模型 对 关键 参数 的 取 值 具有 很 大 随机 性 ， 本 文 为 了 客观 评估 改进 算法 的 性 能 ， 利 用 表 1 所 示 的 两 组 基 
将 利用 OSHHO 算法 对 SVM 的 关键 参数 寻 优 ， 再 利用 优化 后 准 函 数 进 行 寻 优 测试 。 一 组 基准 函数 为 fF1~F4， 为 单 峰 函数 
的 SVM 模型 对 新 冠 肺 炎 疫 情 传 播 趋势 进行 预测 。 类 型 ， 另 一 组 基准 函数 为 FS~F8， 为 多 峰 函 数 类 型 。 基 准 函 
3.2 ”预测 算法 OSHHO-SVM 数 定义 域 对 应 于 个 体位 置 中 每 个 维度 的 搜索 范围 。 选 取 引 言 
利用 OSHHO 算法 对 SVM 的 惩罚 因子 c 和 宽度 参数 o。 ”部 分 中 提 及 的 DEHHO 算法 00、CHHO 算法 03] 和 OBLHHO 
进行 自 适 应 寻 优 ， 构 建 最 优 支 持 向 量 机 预测 模型 ， 实 现 对 新 算法 09 与 本 文 的 OSHHO 算法 进行 性 能 对 比 。 实 验 在 相同 的 
冠 肺 炎 疫 情 传播 的 预测 。 预 测算 法 过 程 如 下 : 硬件 平台 下 进行 ， 实 验 环 境 为 Matlab， 算 法 的 公有 参数 设置 
步骤 1: 确定 支持 向 量 机 SVM 的 输入 量 以 及 惩罚 因子 c ”为 :种 群 规模 为 30, 位 置 维度 为 100, 算法 迭代 总 次 数 为 500， 
和 核 函 数 宽度 参数 o 的 搜索 范围 ， 确 定 OSHHO 算法 的 初始 算法 运行 次 数 为 20 次 。 
参数 : 种 群 规模 W、 最 大 迭代 次 数 Tmax。 表 2 是 四 种 算法 在 八 个 基准 函数 上 的 寻 优 均值 和 标准 方 


步骤 2: 按 比 例 将 原始 数据 集 划 分 为 让 
本 ， 并 以 滑动 窗 
步骤 3: 根据 SVM 模型 的 结构 对 哈 旦 
编码 。 以 均 方 误差 函数 MSE 定义 评估 个 


函数 ， 有 具体 为 


其 中 ,5 为 样 
实际 确诊 人 数 。 


本 量 ， yi 为 SVM 模型 


了 


口 方法 对 训练 集 数据 进行 重 


练 村 
构 。 


本科 


[测试 样 


体位 


>o- DA 


Fitness = 
ss 


斯 鹰 个 体位 置 进行 
劣 的 适应 度 


置 优 
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差 结 果 。 对 于 单 峰 基 准 函数 Fl1~F4， 由 于 高 


了 本 文 的 OS 


高 乡 


HHO 算法 在 F2 上 可 以 找到 最 优 


C7) 


4 的 确诊 人 数 预 测 值 ，yj 为 


算法 几乎 无 法 获 
算法 虽然 未 在 环 
高 的 。F3 是 一 


种 病态 复杂 的 


而 F4 是 一 种 
精度 上 讲 ， 


返 代 后 


OSHHO 算法 
量 级 ， 说 明 算法 依然 具备 竞争 
期 进行 无 c 点 变异 之 后 , 可 以 跳 离 阶梯 状 的 局 部 最 优 处 ， 


三 


取 最 优 解 ， 都 售 
1 上 得 到 最 优 


留 在 


解 ， 


局 部 最 优 


与 四 


但 由 于 阶梯 较 多 ， 在 寻 优 末期 


其 他 三 种 


解 处 。OSHHO 
解 ， 但 其 收敛 精度 是 算法 中 最 
单 峰 函数 ， 多 数 函 数 表 现 较 差 ; 
有 阶梯 状 的 特殊 函数 , 极 易 陷入 局 
种 算法 中 得 到 最 人 
力 。 这 也 证 实 帮 


部 极 值 点 。 从 
CE 解 仅 相 差 三 个 
FE OSHHO 算法 的 


算法 即 停止 了 寻 优 进化 。 
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表 1 基准 函数 说 明 论 最 优 值 , 但 其 收敛 精度 是 所 有 算法 中 最 高 的 。 尤其 对 于 F8 
Tab.1 Reference function description 函数 ， 其 分 布 特征 是 存在 多 个 波峰 极 值 ， 同时 距离 理论 最 [ 
基准 函数 函数 名 称 搜索 区 域 理论 最 优 解 解 较 远 的 局 部 极 值 点 ， 本 文 的 OSHHO 算法 依然 可 以 通过 无 
Fl Sphere [-100,100 0 迹 o 变异 扩大 种 群 搜 索 范围 ， 最 大 限度 地 跳 离 局 部 极 值 点 ， 
F2 Quartic [-1.28,1.28] 0 靠近 理论 最 优 解 ， 较 好 地 平衡 了 搜索 与 开发 间 的 关系 。 而 标 
F3 Rosenbrock [-30,30] 0 准 方差 值 越 小 ， 证 明 OSHHO 算法 在 不 同 基准 函数 上 寻 优 具 
F4 Step [-100,100 0 有 更 好 的 稳定 性 ， 在 具有 不 同 波峰 形态 和 局 部 极 值 的 函数 上 
F5 Ackley [-32,32] 0 寻 优 具有 更 好 的 适应 性 。 
F6 Rastrigin [-5.12,5.12] 0 图 5 是 算法 的 收敛 曲线 图 。 从 500 次 迭代 的 收敛 精 
F7 Griewank [-600,600 0 度 上 看 ,OSHHO 算法 的 收敛 性 能 要 优 于 三 种 对 比 算法 。 
F8 Schwefel [-10,10] 0 此 外 ， 对 于 F1、F5、F6、F7、F8 五 个 基准 函数 而 言 ， 
表 2 算法 在 基准 函数 上 寻 优 的 均值 结果 OSHHO 算法 的 收敛 曲线 具有 较 明 显 的 转折 坠落 ， 这 种 
对 于 多 峰 基 准 函 数 FS~F8， 此 类 函数 通常 具有 多 个 局 部 。 转折 在 不 同 函数 上 所 处 的 迭代 次 数 各 不 相同 ， 这 说 明 算 
极 值 点 ， 其 位 置 分 布 也 各 不 相同 ， 对 于 寻 优 算法 要 求 更 高 ， 法 的 寻 优 精度 在 大 幅 提 升 ， 得 益 于 OSHHO 算法 的 无 迹 
极 易 落 入 局 部 极 值 。 根据 结果 ，OSHHO 算法 在 F6、F7 两 个 o 变异 机 制 能 够 帮助 种 群 个 体 拓展 解 的 搜索 空间 范围 ， 
基准 函数 上 得 到 了 理论 最 优 值 ， 在 FS、F8 上 虽然 未 得 到 理 。 得 到 更 多 的 可 行 解 。 


Tab.2 The algorithm optimizes the mean value result on the benchmark function 


均值 
算法 
Fl F2 F3 F4 F5 F6 F7 F8 
DEHHO 4.12E-76 3.09E-69 2.23E-89 8.45E-67 6.31E-79 1.13E-88 5.85E-101 4.83E-97 
CHHO 1.12E-87 6.35E-75 7.45E-102 4.61E-79 5.11E-92 7.15E-97 5.59E-121 5.87E-119 
OBLHHO 6.89E-101 3.98E-111 4.12E-119 3.15E-102 4.63E-111 5.98E-132 5.64E-141 4.49E-138 
OSHHO 4.28E-113 0 5.37E-112 6.03E-99 7.37E-132 0 0 5.37E-142 
算法 
Fl F2 F3 F4 F5 F6 F7 F8 
DEHHO 6.30E-76 2.76E-69 4.25E-89 4.45E-67 3.81E-79 5.63E-88 2.45E-101 4.38E-97 
CHHO 4.32E-87 1.18E-75 4.75E-102 6.34E-79 5.02E-92 5.35E-97 3.19E-121 7.47E-119 
OBLHHO 5.19E-101 5.68E-111 5.08E-108 5.11E-102 4.11E-111 6.88E-132 4.04E-141 5.58E-138 
OSHHO 0 0 7.38E-112 6.78E-99 3.31E-132 2.04E-149 0 4.18E-144 
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图 5 算法 的 收敛 曲线 


Fig.5 Convergence curve of algorithm 
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录用 定稿 

为 了 进一步 分 析 
算法 全 局 搜索 能 力 的 影响 ， 
以 图 示 方 式 展示 最 优 解 个 


陈 强 ,， 等 : 


到 6 所 示 。 该 图 以 基准 


函数 Rosenbroc 


行 测试 ， 总 共 适 
其 他 基准 函数 上 进行 。 根 


搜索 区 域 广泛 ， 全 局 搜索 
解 的 区 域 广泛 寻 优 ， 保 持 


局 


尺 200 次 ， 相 关 实 验 也 可 以 扩展 到 
据 测 试 结 果 ， 搜 索 个 体 在 欠 代 初期 
能 力 较 强 ， 能 够 充分 的 在 远离 候选 


伪 对 立 学 习 和 伪 反 射 学 习 机 制 对 HHO 
以 单个 种 群 个 体 的 搜索 轨迹 为 例 ， 
体 与 搜索 个 体 的 搜索 轨迹 对 比 图 ， 


k 和 Ackley 搜索 的 


着 较 好 的 多 样 性 ， 而 


随 着 迭代 的 进 


行 ， 后 期 搜索 个 体 更 加 倾 


向 于 


局 部 开采 和 精 儿 


搜索 ， 并 不 断 


逼近 候选 解 ， 提 高 了 收敛 


速度 和 寻 优 精度 。 


通 
在 HHO 算法 的 软 包 围 和 
反射 学 习 机 制 对 搜索 和 开 


硬 包 围 中 分 别 引 


以 上 结果 表明 ， 
| 入 伪 对 立 学 习 和 僵 


发 过 程 的 均衡 优 
法 的 全 局 寻 优 能 力 和 效率 。 


化 后 ， 能 够 提升 算 


最 优 解 个 体 的 寻 优 曲 线 CS 
4 一 最 优 解 个 体 的 寻 优 曲线 
0 J 休 0 3 由 
| 
、 、 
|| 本 ,| | SA 
性 | 搜索 个 体 的 寻 优 曲线 
-| 搜索 个 体 的 寻 优 曲线 -| 
0 20 40 60 80 100 120 140 160 180 200 本 0 20 40 60 100 120 
连 代 选 代 
(a)Rosenbrock 函数 (b)Ackley 函数 
图 6 个 体 搜索 轨迹 


Fig.6 Individual search trajectory 


测试 函 
如 图 7 是 四 种 算法 在 种 群 


400 次 时 种 群 个 体 在 搜索 空 


多 样 性 。DEHHO 和 CHHO 算法 


始终 比较 分 散 ， 尤 其 进入 


数 F4 在 不 同 迭 代 次 数 时 种 群 个 体 


的 分 布 情况 。 


规模 为 30， 和 迭代 100 次 、200 次 和 


E 间 内 的 分 布 。 


碗 代 后 


极 值 点 而 无 法 脱离 ， 离 最 


改进 策略 还 是 无 法 根本 改 


基部 分 个 体 已 聚集 于 最 优 


法 的 寻 优 精度 。 本 文 的 O 
种 群 进化 没有 过 早 聚 集 在 
索 与 开发 间 的 均衡 以 及 跳 
略 的 混合 改进 还 是 起 到 
经 逼近 于 最 优 解 附近 。 


进一步 观察 在 固定 目 
法 达到 该 寻 优 精度 所 需 的 


判定 标准 为 


两 种 算法 , OBLHHO 算法 全 局 收敛 性 
解 区 域 ， 对 立 学 习 和 正 交 学 习 使 得 
HHO 算法 中 种 群 能 够 更 好 的 实现 信息 交互 ， 


在 三 个 迭代 时 期 的 种 群 个 体 
期 ， 部 分 个 体 逐 步 陷入 局 部 
优 解 区 域 依然 较 远 ， 


个 体 分 布 反映 种 群 


这 表明 单一 的 


善 HHO 算法 


的 寻 精度 。 


相 比 这 


尤其 在 迭代 后 


很 好 地 提升 了 算 


SHHO 算法 则 


具有 最 好 的 聚集 性 ， 


局 部 最 优 区 域 ， 


对 种 群 多 样 性 、 搜 


离 局 部 极 值 方 


田 ， 该 


和 法 通过 多 策 


标 基准 


[FEFl/F <es,F, #0 
IE-Eles,F,=0 


其 中 : F4 为 实际 求解 最 优 解 ， 


Fr 为 基 ? 


信函 数理 论 最 优 解 。 


预期 效果 ， 和 迭代 后 期 ， 多 数 个 体 已 


函数 的 收敛 精度 时 ， 不 同 算 
迭代 次 数 和 相应 的 寻 优 成 功率 。 寻 
成 功率 为 寻 优 成 功 次 数 与 实验 次 数 的 比值 ， 


而 寻 优 成 功 的 


(28) 


设置 函数 维度 为 20， 最 大 迭代 数 为 500， 表 3 是 算法 的 


测试 结果 。 当 从 代 500 次 
代 次 数 填 为 500。 
500 次 均 未 达到 目标 寻 优 


仍 未 到 达 


DEHHO 算法 在 函数 F2、F3、F5 


精度 ， 所 以 其 在 这 


寻 优 成 功率 都 是 0， 算 法 的 寻 优 性 能 还 


OSHHO 算法 在 fF2、F6、 
同时 利 
和 随机 无 迹 c 变异 机 制 更 
了 种 群 向 着 最 优 解 区 域 的 


F7 上 均 达 到 


所 需 寻 优 精度 时 ， 则 将 迭 
上 从 代 了 
个 基准 函数 上 


有 待 提升 。 本 文 的 


100% 的 寻 优 成 功率 ， 
了 最 少 的 迭代 次 数 ， 说 明 改 进 算法 


好 的 引导 了 种 群 


所 采用 的 伪 学 习 
的 寻 优 方向 ， 加 快 


转移 速度 。 


本 文 提 出 的 算法 不 仅 寻 


优 成 功率 更 加 稳定 ， 同 时 
民 次 数 也 是 最 少 的 ， 说 明 


其 得 到 


标 寻 优 精度 时 所 需要 的 和 迭 


算法 处 理 复 杂 多 峰 问题 时 ， 也 能 够 


效 提升 寻 优 精度 和 收敛 


速度 。 


搜索 空 
算法 最 大 友 代 次 数 为 Tnax， 空 间 维 
的 时 间 复 杂 度 为 OCWWxq)。 根 据 随 机 无 迹 c 点 变异 机 制 ， 
OSHHO 算法 需要 产生 24+1 个 随机 点 并 对 其 进行 排序 ， 


1)OSHHO 算法 的 


ChinaXiv 合 作 期 刊 


一 种 基于 随机 无 迹 o 变异 的 改进 HHO 算法 及 应 用 


时 间 复 杂 度 分 析 


第 39 卷 第 9 期 


群体 智能 优化 算法 的 时 间 复 杂 度 主要 取决 于 种 群 规模 、 


s 间 维 


杂 度 为 O((mitn2)xqdxTia)。 


速 排序 对 随机 点 进行 排序 ， 可 知 该 
O((Q24d+1)+(24+1)lg(24+1)), 
的 种 群 个 体 数 量 分 别 为 n! 和 n2， 上 是 


度 以 及 算法 的 最 大 迭代 次 数 。 令 种 群 规模 为 N， 
度 为 4。 算 法 初始 化 过 程 


以 快 


阶段 的 时 间 复 杂 度 为 


令 进 行 伪 对 立 学 习 和 
n1、n2<N， 内 
综 上 ，OSHHO 算法 的 最 差 复杂 度 为 
OVxdiHO(CCdHDHCdHDlegCdaHD)HO(UnHDp)xdxThao=OGCVxdx7hao ， 


这 表明 改进 算法 并 未 在 HHO 算法 的 基 


擅 反 射 学 习 
该 过 程 的 时 间 复 


基础 上 增加 时 间 代 价 ， 


与 原 算法 HHO 是 保持 一 致 的 。 


0 = 


1 
选 代 100 次 


-i 
选 代 200 次 


加 
和 迭代 400 次 


(a)DEHHO 算法 
100 io 100 
| 本 
1 00 100 
选 代 100 次 迁 代 200 次 选 代 400 次 
(b)CHHO 算法 
100 100 100 
00 -0 ai 
办 代 100 次 迁 代 200 次 和 迭代 400 次 
(c)OBLHHO 算法 
io0 
i i 
1 5 而 
迁 代 100 次 迁 代 200 次 迁 代 400 次 
(d)OSHHO 


区 | 


7 和 迭代 中 的 种 群 个 体 分 布 


Fig.7 。 Population individual distribution in iteration 
表 3 ， 寻 优 成 功率 及 收敛 迭代 数 


Tab.3 Optimization success rate and convergence iteration 


管 法 寻 优 成 功率 /2% 
Fl F2 F3 F4 FS F6 F7 F8 
DEHHO 93.4 0 0 78.9 0 92.1 90.5 93.2 
CHHO 91.1 86.3 96.7 97.4 98.4 97.9 94.9 94.2 
OBLHHO 92.4 95.7 96.2 98.1 99.2 98.8 95.3 94.0 
OSHHO 99.5 100 99.6 99.1 99.3 100 100 98.9 
管 法 寻 优 迭代 次 数 /次 
Fl F2 F3 F4 FS F6 F7 F8 
DEHHO 498 500 500 489 500 478 469 493 
CHHO 409 401 404 435 489 398 421 445 
OBLHHO 487 389 399 338 410 376 367 401 
OSHHO 328 266 231 208 395 319 332 389 
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2)OSHHO 算法 的 收敛 性 分 析 

定理 1 OSHHO 算法 中 种 群 个 体位 置 更 新 解 集 合 
fxx(D, 它 0} 是 满足 在 有 限 状 态 空间 内 的 齐 次 Markov 链 。 

证 明 OSHHO 算法 的 迭代 寻 优 过 程 可 知 ， 每 次 迭代 
中 最 优 解 的 更 新 均 局 限于 有 限 状 态 空 间 中 。 针 对 解 更 新 所 采 
用 的 改进 策略 伪 对 立 和 伪 反 射 学 习 策 略 、 能 量 因子 非 线性 周 
期 性 调整 以 及 随机 无 迹 o 点 变异 机 制 都 满足 相互 独立 条 件 ， 
即 针对 哈里 斯 诬 个 体 寻 优 位 置 最 新 中 是 一 种 逐步 更 新 机 制 


实际 确诊 人 数 / 个 


当 次 欠 代 (0 时 的 种 群 仅 接收 前 代 友 代 (zcD 时 种 群 的 状态 信息 
以 触发 更 新 过 程 , 因此 , 算法 的 解 集合 fxxv(D, 之 0} 是 满足 在 有 样本 天 数 /天 
限 状态 空间 内 的 齐 次 Markov 链 。 图 8 日 增 确诊 数据 
引 理 1 OSHHO 算法 的 状态 空间 内 ， 存 在 使 条 件 pij<x Fig.8 Daily diagnostic data 
成 立 的 常量 jE(0,1),pij 为 状态 空间 内 状态 i 转移 至 状态 j 的 表 4 训练 集 数据 的 重 构 
概率 。 Tab.4 Reconstruction oftraining set data 
定理 2 OSHHO 算法 收敛 到 全 局 最 优 的 概率 为 1 。 ] 维 输入 向 量 一 维 输入 向 量 
证 明 给 定 搜索 精度 s>0, 则 Markov 链 的 序列 状态 空间 x1,x2,x3,x4 x5 
可 划分 为 X2,X3,X4,XS x6 
Q={d| f(D)- fl<e} Xx3,X4,x5,x6 和 
2 ={d7OD- 记 3} 本 ... 
oe = min{ f (x;(D),i=1,2,...,N} X71,X72,xX73,x74 X75 
种 群 从 1 友 代 至 上 代 时 , 令 无 法 到 达 搜 索 精度 的 概率 选择 标准 支持 向 量 机 SVM 模型 、 标 准 哈里 斯 雇 算 法 优 
为 2， 则 根据 引 理 1 有 : 化 SVM 模型 HHO-SVMDGI、 粒 子 群 优化 神经 网 络 PSO-BPD2] 
pe 以 及 基于 OBLHHO 算法 Ms 优化 SVM 模型 OBLHHO-SVM 
三 预测 模型 进行 实验 对 比分 析 。 与 标准 SVM 模型 对 比 可 以 验 
再 根据 波 莱 尔 - 坎 泰利 引 理 ， 可 得 证 OSHHO 算法 优化 SVM 的 可 行 性 ， 与 HHO-SVM 、 
ple FOOD- 所 习 =0 OBLHHO-SVM 模型 对 比 可 以 验证 OSHHO 算法 的 优越 性 ， 
已 知 结论 : 对 于 随机 序列 60， 若 存在 09， 使 得 任意 仿 0 与 PSO-BP 模型 对 比 则 可 以 验证 所 选 支持 向 量 机 模型 对 于 新 
时 ， 满 足 : 冠 疫 情 传播 预测 是 否 有 效 有 行 。 PSO-BP 预测 模型 中 , 学 习 因 
PM Bi)=0 子 设置 为 2， 惯性 权重 最 小 值 为 0.1， 最 大 值 为 0.9， 神 经 网 
则 表明 随机 序列 {9()} 能 收敛 于 0。 络 结构 为 10-6-1， 即 输入 层 神 经 元 数量 为 10， 隐 含 层 神经 元 
综 上 结论 可 知 ，OSHHO 算法 收敛 到 全 局 最 优 的 概率 为 1。 ” 数 为 6， 输 出 层 神经 元 数 为 1( 对 应 预测 值 )， 神 经 网 络 训练 
4.2 实验 2: 基于 OSHHO-SVM 模型 的 新 冠 疫情 预测 标 误差 为 1x105， 训 练 次 数 为 1000， 学 习 率 为 0.05 。 
TD) 实验 环境 与 样本 说 明 2) 实 验 结果 
选取 某 地 区 在 2020 年 3 月 10 日 ~6 月 17 日 之 间 100 天 如 图 9 是 五 种 预测 模型 在 测试 样本 上 得 到 的 确诊 人 数 的 
的 日 增 肺炎 确诊 人 数 作为 样本 数据 , 测试 OSHHO-SVM 模型 ”预测 值 与 实际 值 间 的 对 比 情况 ， 不 同 预测 模型 的 曲线 表示 为 
对 于 预测 肺炎 传播 的 性 能 。 该 样本 数据 来 源 于 国家 卫 健 委 虚线， 确诊 人 数 的 实际 值 表示 为 实 线 。 总 体 来 看 ， 将 智能 优 
(http:/wwwnhc.gov.cn/xcs/yqtb/list_gzbd.shtml) 每 日 疫情 通报 数据 ”化 算法 与 支持 向 量 机 模型 或 BP 神经 网 络 相 结合 ， 对 预测 模 


中 的 日 增 确诊 人 数 。 图 8 是 该 段 时 间 区 间 内 经 预 处 理 后 得 到 ” 型 的 准确 率 有 一 定 程度 的 提升 ， 在 预测 新 冠 肺 炎 疫 情 传播 问 
的 每 日 确诊 人 数 变化 情况 。 从 2020 年 3 月 10 日 之 后 , 疫情 。 题 上 具有 可 行 性 。 从 预测 值 与 实际 值 曲线 的 拟 合 程度 上 看 ， 
传播 逐渐 处 于 爆发 期 ， 日 增 人 数 迅 猛 增加 。 虽 然 在 政府 防 控 ” 模型 具有 一 定 差异 性 。 标准 SVM 模型 、 PSO-BP 模型 、HHO- 
介入 严禁 人 口 流 动 和 人 人 佩戴 口罩 等 有 效 防 疫 措施 实施 后 ， SVM 模型 以 及 OBLHHO-SVM 模型 在 预测 值 的 拟 合 度 上 都 
后 期 日 增 人 数 相 应 降低 ， 但 依然 还 存在 局 部 性 爆发 趋势 。 从 不同 程度 偏差 ， 利 用 OBHHO 算法 优化 的 支持 向 量 机 预测 
整个 曲线 趋势 可 以 看 出 ， 确 诊 人 数 变 化 是 典型 的 非 线 性 、 非 ” 模型 在 训练 后 ， 预 测 值 与 实际 值 具 有 最 大 程度 的 拟 合 ， 大 部 
稳定 的 分 布 特征 ， 局 部 位 置 显现 多 个 峰值 。 传 统 SVM 模型 ” 分 日 期 中 的 预测 准确 率 达 到 100%, 说 明 算法 对 SVM 的 惩罚 
对 这 类 数据 存在 拟 合 困难 ， 预 测 准确 率 低 ， 需 要 对 SVM 的 ” 参数 和 核 函 数 宽度 参数 进行 调整 优化 后 ， 算 法 具有 不 断 向 最 
关键 参数 : 惩罚 因子 c 和 宽度 参数 o 进行 迭代 优化 ， 才 能 使 优 解 逼近 的 能 力 ， 算 法 的 适应 度 在 不 断 提 升 。 标 准 SVM 模 
预测 模型 具有 更 高 的 准确 率 。 型 的 预测 偏差 相对 较 大 , 模型 泛 化 能 力 和 稳定 性 都 有 待 提升 。 

将 样本 数据 以 3:1 的 比例 划分 为 训练 集 和 测试 集 ， 模 型 HHO-SVM 模型 和 OBLHHO-SVM 模型 在 部 分 日 期 上 确诊 人 
独立 运行 20 次 ， 样 本 前 75 天 数据 为 训练 集 ， 后 25 天 数据 ” 数 预 测 值 与 实际 值 极为 接近 ， 拟 合 度 优 于 PSO-BP 模型 ， 这 
为 测试 集 。 相 关 参 数 中 ,种群 规模 设置 N=50, 算法 最 大 迭代  ” 表明 算法 的 寻 优 精度 还 有 提高 空间 ,综合 比较 , OSHHO-SVM 


三 


次 数 Tmax=400。SVM 模型 中 ,惩罚 因子 c 和 宽度 参数 c 的 搜 。 利用 改进 的 哈里 斯 座 优 化 算法 提高 了 模型 的 预测 精度 和 泛 化 
索 范 围 设置 为 [2-5,29]。 能 力 ， 具 有 性 能 优势 。 
于 新 冠 疫 情 的 爆发 具有 连续 性 和 随机 性 ， 若 预测 时 间 选择 平均 绝对 误差 M4E、 平均 相 对 误差 RMSE 和 拟 合 度 


间隔 过 小 ， 会 导致 模型 的 学 习 程度 不 足 ， 而 预测 时 间 间隔 过 ”因子 R? 三 个 指标 评估 算法 的 预测 性 能 。 相 关 定义 如 下 

大 ， 又 会 影响 到 模型 对 邻近 时 段 内 疫情 传播 的 学 习 。 故 针对 ji 天 
前 75 天 的 训练 集 数 据 和 后 25 天 的 测试 集 数据 ， 利 用 滑动 窗 5 

口 法 对 训练 集 数据 重 构 ， 窗 口 设置 为 5， 且 每 次 滑动 单元 数 
设 为 1。 具 体重 构 结果 如 表 4 所 示 。 


(30) 
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全 局 搜索 能 力 进行 了 改进 。 结 合 2019 年 爆发 的 新 冠 肺炎 疫 


情 传播 情况 ， 利 用 OSHHO 算法 对 支持 向 量 机 模型 SVM 进 
行 优化 ， 建 立 了 面向 新 冠 肺 炎 疫 情 传播 的 预测 模型 OSHHO- 
SVM。 并 利用 若干 样本 数据 证 明 预 测 模型 得 到 的 拟 合 度 、 预 


六 


其 中 , 7 为 选择 的 样本 数量 ,六 为 样本 实际 值 ， 为 水 肥 用 量 
预测 值 ，yave 为 实际 值 均值 ，yiave’ 为 预测 值 均值 。 评 估 指 标 
中 , MAE 用 于 评估 预测 值 误 差 的 实际 情况 ，MA4E 越 小 , 表明 
算法 精度 越 高 ; RMSE 用 于 评估 预测 可 信和 度 ,， RMSE 越 小 , 表 
明 算 法 精度 更 高 ， 结 果 更 可 信 ; 拟 合 度 因子 则 评估 了 预测 值 


与 实际 值 之 间 的 拟 合 度 ，R?<1， 值 越 大 ， 拟 合 越 高 。 
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图 9 预测 值 与 实际 值 对 
Fig.9 Comparison between predicted value and actual value 
如 表 5 是 以 上 三 个 评价 指标 的 对 比 结果 。 对 于 M4E 和 

RMSE，OSHHO-SVM 模型 得 到 的 值 是 所 有 模型 中 最 小 的 ， 

正明 改进 算法 的 预测 精度 是 最 高 的 。 对 于 及， 标准 SVM 模 


ENY 


iv 


YH 


型 离 1 相差 最 远 ,说 明 模型 预测 值 与 实际 值 的 拟 合 程度 较 差 。 


= 
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0 种 启发 式 优化 算法 结合 支持 向 量 机 的 预测 模型 都 在 不 同 程 
度 上 提高 了 数据 拟 合 度 ， 即 相 应 提升 了 R? 值 , 证 明 通 过 相应 
的 优化 机 制 针 对 惩罚 参数 和 核 函 数 宽度 参数 的 优化 ， 能 够 提 
高 SVM 的 预测 能 力 。OSHHO-SVM 模型 在 R2 值 上 与 极限 值 
1 仅 相 差 0.0184, 表明 改进 策略 在 提升 SVM 的 模型 泛 化 能 
以 及 预测 精度 上 具有 现实 可 行 性 。 


IT 


表 5 评估 指标 结 
Tab.5 _ Evaluation index results 
预测 模型 评估 指标 
MAE RMSE R2 
标准 SVM 模型 1.5745 3.7373 0.8416 
PSO-BP 模型 1.5632 3.4092 0.8509 
HHO-SVM 模型 1.5028 3.1378 0.9023 
OBLHHO-SVM 模型 1.3113 2.7093 0.9668 
OSHHO-SVM 模型 1.1937 2.6736 0.9816 
如 图 10 是 模型 收敛 性 的 对 比 结果 。 依 据 适 应 度 函 数 的 


定义 ， 适 应 度 值 越 小 ， 表 明 预 测 模型 求解 的 确诊 人 数 误差 越 
小 。 从 整个 收敛 曲线 趋势 看 ，OSHHO-SVM 模型 得 到 最 佳 的 
适应 度 值 ,表明 最 终 的 预测 精度 是 最 高 的 。 而 在 收敛 速度 上 ， 
OSHHO-SVM 模型 在 58 次 迭代 后 得 到 最 佳 适 应 度 ， 表 明 此 
时 通过 OSHHO 算法 求解 到 了 最 优 的 惩罚 参数 和 核 函数 宽度 
参数 , 使 SVM 具备 了 最 优 的 预测 性 能 。 比较 而 言 , OBLHHO- 
SVM、HHO-SVM 和 PSO-BP 三 种 预测 模型 分 别 在 经 过 88、 
95、101 次 迭代 后 收敛 于 最 佳 适应 度 值 ，OSHHO-SVM 相对 
这 三 种 模型 在 收敛 速度 上 分 别提 高 了 34.09%、38.94%、 
42.57%。 


5 ”结束语 

本 文 提出 了 一 种 融合 无 迹 变异 的 改进 哈里 斯 鹰 优化 算法 
OSHHO。 改进 算法 分 别 利 用 伪 对 立 和 伪 反 射 学 习 机 制 、 能 量 
因子 非 线性 周期 性 调整 机 制 以 及 随机 无 迹 o 点 变异 机 制 对 


HHO 算法 的 寻 优 效率 、 全 局 搜索 与 局 部 开发 能 力 的 均衡 以 及 


测 精 度 都 得 到 了 稳步 提升 。 进 一 步 的 研究 可 集中 在 种 群 划分 
问题 上 ,针对 不 同 优 劣 的 种 群 个 体 (如 优 、 中 、 差 ) 利 用 不 同 机 
制 的 位 置 进化 策略 ， 加 快 算法 围绕 最 优 解 的 渐 近 模式 ， 提 升 
算法 整体 的 寻 优 效率 。 
1 一 一 PSO-BP 模 型 
s =- = = "HHO-SVM 模 型 
0.75 一 -= “OBLHHO-SVM 模 型 
一 一 = “OSHHO-SVM 模 型 
其 
全 ”0.5 
ba) ~ 
025 | \ Pe 
Me 
0 1 1 
0 100 200 300 400 
迭代 /次 
图 10 模型 收敛 性 比较 
Fig. 10 Comparative model convergence 
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